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Abstract: Mounds constructed by plateau zokors, which is widely distributed in alpine meadows 
significantly modified plant community structure. However, the variations of plant community structure 
under the disturbance of plateau zokor-made mound are less concerned. Therefore, we investigated the 
responses of plant community on zokor-made mound of different years (1 a and 3-4 a), and compared 
with undisturbed sites (no mound) in an alpine meadow in the eastern Qinghai-Tibetan Plateau (QTP), 
China. Species richness, coverage and Simpson diversity index were all significantly reduced by the 
presence of zokor-made mound, but plant heights were significantly increased, particularly in grasses and 
sedges. Several perennial forage species showed an increased importance value and niche breadth, 
including Koeleria macrantha, Elymus nutans and Poa pratensis. The effect of zokor-made mound on niche 
overlap showed that more intense interspecific competition produced a greater utilization of 
environmental resources. And this interspecific niche overlap was strengthened as succession progressed. 
The bare mound created by zokor burrowing activities provided a colonizing opportunity for 
non-dominant forage species, resulting in abundant plant species and plant diversity during the succession 
period. We concluded that presence of zokor-made mound was conducive to regeneration and vitality of 
plant community in alpine meadows, thus improving their resilience to anthropogenic stress. 
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1 Introduction 


The Qinghai-Tibetan Plateau (QTP), China is the world's largest alpine pastoral ecosystem. 
It provides ecological benefits of carbon sequestration, soil and water resource protection, 
aesthetic recreation and tourism, and economic benefits of pastoral production (Feng et al., 
2010; Wang et al., 2016; Hopping et al., 2018; Zhang et al., 2019; Dong et al., 2020; Liu et al., 
2020). However, these meadows are threatened by increased livestock grazing and human 
population growth (Harris, 2010; Nyima, 2017) and a changing climate that is thawing 
permafrost and facilitating 
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population growth of burrowing rodents (Xue et al., 2009; Harris, 2010; Xue et al., 2017), 
resulting in varying degrees of degradation. Alpine meadow degradation threatens pastoral 
livelihoods and ecological function via soil erosion and nutrient loss (Li et al., 2013; Lu et al., 
2017; Liu et al., 2018; Wang et al., 2020; Yuan et al., 2020). 

The plateau zokor (Myospalax baileyi) is a subterranean rodent endemic to the alpine meadow 
of the QTP, China (Wang et al., 2008; Hu et al., 2017). This species is often described as a major 
pest due to its foraging, digging and mound-making activities producing highly visible changes to 
the meadow landscape. Its impact on plant diversity, biomass, soil organic carbon, temperature 
and soil water holding capacity is usually characterized as negative (Wang et al., 2008; Li et al., 
2009; Zhang et al., 2014; Hu et al., 2017; Chu et al., 2020). Enzyme activity and microbial 
biomass are lower in the zokor-made mound than in the surrounding undisturbed soil (Li et al., 
2009). Nevertheless, recent studies have suggested the positive ecological benefits of the plateau 
Zokor, including increased soil water availability, soil organic carbon, and diverse plant 
communities across the meadow (Wang et al., 2008; Wu et al., 2017). Seed germination and 
seedling survival are both higher on the zokor-made mound (Bao et al., 2016; Hu et al., 2017). 
Moderate habitat disturbance by plateau zokor may stimulate species recruitment and increase 
meadow plant diversity, and therefore improve regeneration after a meadow has been overgrazed 
(Wu et al., 2015; Hu et al., 2017). 

There have been few studies on niche characteristics of plant communities in alpine meadows 
disturbed by plateau zokor (Hu et al., 2017). To address this gap, we studied plant community 
structure on zokor-made mound in different ages. Our aim was to (1) characterize the plant 
community succession after disturbance by mound construction in meadow; (2) determine the 
responses of niche characteristics of plant community to mound construction; and (3) compare the 
effects of zokor-made mound on the recovery of plant diversity. 


2 Materials and methods 


2.1 Study area 


The research area was located at the Hongyuan County Alpine Meadow Ecosystem Research 
Station of the Southwest Minzu University, China. Hongyuan County lies on the eastern margin 
of the QTP (31?50'-33?22'N, 101?51'-103?23'E). Annual precipitation is 650—730 mm falling 
mostly between May and August. Average altitude is 3500 m a.s.l., annual mean temperature is 
1.1?C and mean annual duration of snow cover is 76 d (Li et al., 2011). 


2.2 Sampling and measurements 


Plant coverage, height and frequency were recorded for each plant species in nine replicate plots 
of three mound treatments in September 2011. Treatments included the control of undisturbed 
meadow (no mound), and zokor-made mound aged one year (1 a) and three to four years (3-4 a) 
since construction, with plots selected randomly from the meadow. Mound had been abandoned 
by zokors at the time of study. The 1 a mound was constructed in August 2010, and at time of 
study had a surface area of 0.6 mí and a height of 80-100 mm. The 3-4 a mound was constructed 
on or before 2008 and at time of study had a surface area of 0.4 m? and a height of 30-50 mm. 
Control plot of 1.0 m? was selected on site that showed no visible signs of mound construction. 
Plant species were classified into four functional groups for further analysis: grass, sedge, legume 
and forb (Hu et al., 2017). 


2.3 Data analysis 

Importance value (IV) of each species was calculated by the following equation: 
IV=(C,+-A;+F,)/3, (1) 

where C,, H, and F, are the relative coverage, relative height and relative frequency, respectively 


(Hu et al., 2017). 
Three diversity indices were calculated: Shannon—Wiener diversity index (H), Simpson 
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diversity index (D; Lande, 1996) and Pielou evenness index (E; Pielou, 1966). The formula of 
calculation of each index is given below: 


rs om Q) 
- 2 

DW 3) 

E=H/InS, (4) 


where P; is the importance value; and S is the plant species richness (Fang et al., 2009). 

Six indices were used to assess the niche characteristics of plant community: (1) Levins niche 
breadth index (NBr; Hurlbert, 1978); (2) Shannon—Wiener niche breadth index (NBsw; Colwell 
and Futuyma, 1971); (3) Hurlbert niche breadth index (NBy; Hurlbert, 1978); (4) Levins 
coefficient of niche overlap (O; Hurlbert, 1978); (5) ecological attributes (^O;) and (6) 
ecological response rate (R; Zhang, 2018). The equations of each index are given below: 


NB, = 1/ (ny P) (5) 
pat 
NB,, = -2 BPR; , (6) 
Ja 

NBg-(nNB;-1)/(n-1), (7) 
O, = ‘os BE / $3 PI (8) 

j=l j=l 
AO, = O, O; = 2; O) , (9) 

k=1 i=l 
r R-NB/AOg; (i-k), (10) 


where P;=n;/Ni; N, = » "i; ; nj is the dominance of the i" population in the j" resource state; Pj 
is the proportion of the dominance of the ^ population in the j” resource state to the total 
dominance of the i” population in all resource states; n is the number of resources states; Py is the 
proportion of the dominance of the p population in the j” resource state to the total dominance of 
the k^ population in all resource states; O; is the amount of occupied resources of population k by 
population i; Ox; is the amount of occupied resources of population i by population k. If O> Okri, 
the population i is developmental, otherwise, the population k is declining; and NB is the niche 
breadth of plant species. The total niche overlap value is equal to the niche overlap value of all 
species pairs divided by the total number of species pairs in plant community. 

All data were checked for normality and homogeneity of variance. Normality was tested by the 
Kolmogorov-Smirnov Z method (P<0.05) using a nonparametric test. Homogeneity of variance 
was tested by the Levene lest method (P<0.05). One-way analysis of variance (ANOVA) with 
Bonferroni's multiple range test (P<0.05) was performed to assess the effect of different types of 
zokor-made mound (1 a, 3-4 a and no mound) on plant species richness, height, H, E and plant 
functional groups, which met the normality and homogeneity of variance. Data that failed the 
normality and homogeneity assumptions for parametric testing were analyzed by the 
nonparametric Kruskal-Wallis H test (P<0.05). Oj, and O;; of plant species in the community were 
compared by the Kruskal-Wallis H method (P<0.05) in the nonparametric test. Differences 
obtained at the P«0.05 level were considered significant. Niche breadth (NBsw, NB, and NBy) 
and ecological response rate (R-NBsw, R-NB and R-NBy) of plant species were analyzed by the 
principal component analysis (PCA) in CANOCO 5.0 software (Version 5.0 for Windows; 
TerBraak and Smilauer, 2012) to identify the changes in niche characteristics among plant 
communities. 
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3 Results 


3.1 Characteristics of plant community structure 


Zokor-made mound had a reduced plant species richness, coverage and Simpson diversity index, 
but the Shannon-Wiener diversity and Pielou evenness indices were unaffected. Plant height in 
Zokor-made mound increased significantly (Table 1). 

Species richness of sedge and legume decreased significantly compared with no mound, but 
species richness of grass and forb were unaffected (Fig. 1). Grass coverage was greater in 3-4 a 
than both 1 a and no mound treatments, the latter of which was not significantly different. There 
was no significant difference in sedge coverage, but legume coverage of 1 a mound was 
significantly lower. Forb coverage was the highest in 1 a mound and the lowest in 3-4 a mound. 
Grass heights of both 1 a and 3-4 a mounds were significantly greater than that of no mound. 
Sedge height of 3-4 a mound was greater than those of both 1 a and no mound treatments. 
Legume and forb heights declined in both 1 a and 3-4 a mounds (Fig. 1). 


Table 1 Effect of zokor-made mound on plant community structure 


Mound type Statistical result 
Index 

la 3-4a No mound ForH P 
Species richness 18.4424.50*^ 13.56+6.86º 21.56+6.93° F55473.81 0.037 
Coverage (94) 82.56+3.17º 91.56+4.82º 97.444] .33^ H-20.25 «0.001 
Height (cm) 13.48+0.72* 13.50+1.08º 12.43+0.69° Fy 24=4.70 0.019 
Simpson diversity index 0.958+0.005* 0.941+0.013º 0.967+0.001° H=21.08 <0.001 
Shannon- Wiener diversity index 2.355+0.428 1.906+0.623 2.538+0.518 Fy 24=3.40 0.050 
Pielou evenness index 0.810+0.080 0.756+0.077 0.831+0.084 F53472.05 0.151 


Note: Different lowercase letters within the same row indicate significant differences among different mound types at P«0.05 level. H is 


statistic value of the nonparametric Kruskal-Wallis H test. Mean+SE. 
EH | a mound 3—4 a mound [E No mound 


, 610 P=0.296 7 (b) P-0018 54 (6) E P=0.074 
FA 54 a 4 4d 204 a 
LET é b 16 4 
3 34 ab 124 
E P. a 
g 24 84 
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124 (h) P<0.001 
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Fig. 1 Species richness (a-d), coverage (e-h) and height (1-1) for each plant functional group and zokor-made 
mound type. Different lowercase letters indicate significant differences among different types of zokor-made 
mound at P<0.05 level. 
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3.2 Niche breadth of plant species 


Zokor-made mound significantly affected the niche breadth of the 41 observed plant species 
(Table 2). The three formulae (NBsw, NB, and NBy) produced different values but similar trends 
across species (Fig. 2). Dominant species among treatments differed. Mound of 1 a was 
dominated by grasses (Koeleria macrantha, Elymus nutans and Poa pratensis) and forbs 
(Stellaria infracta, Ranunculus tanguticus, Artemisia annua and Potentilla anserina). Mound of 
3-4 a was also dominated by grasses (P pratensis, E. nutans and K. macrantha) and forbs (S. 
infracta and Lancea tibetica). However, in no mound treatment, forbs (Anemone rivularis, 
Saussurea nigrescens, R. tanguticus, P. anserina, Anaphalis lactea and Gentiana ornate), legumes 
(Tibetia himalaica and Oxytropis ochrocephala), one grass (Festuca rubra) and one sedge (Carex 
kansuensis) were dominant (Table 2). 


Table2 Importance value (IV) and ecological attribute value (^O;) for each plant species under different types 
of zokor-made mound 


No. Species Functional Life la 3—4 a No mound 
group cycle IV AO IV AO; IV AO 
SI Koeleria macrantha (Ledeb.) Schult. Grass P 0.067 11.909 0.090 7.422 0.029 6.962 
S2 Elymus nutans Griseb. Grass P 0.048 9.767 0.076 12.948 0.014 —10.019 
S3 Avena fatua L. Grass A 0.007 —10.809 - - - - 
S4 Poa pratensis L. Grass P 0.043 6.442 0.073 15.020 0.013 —10.020 
S5 Festuca rubra L. Grass P 0.016 —0.683 0.010 —8.317 0.051 8.599 
S6 Festuca ovina L. Grass P - - 0.038 —1.328 0.011 —10.134 
S7 Kobresia capillifolia (Decne.) C. B. Sede P 0.011 7.626 ] i 0.044 2442 
Clarke 
gg. dc distigmaricum CSHE] Sedge P — 0.080 0.733 0.054  -4205 0.031 0.686 
T. V. Egorova 
S9 Carex kansuensis Nelmes Sedge P 0.026 1.543 0.087 —0.848 . 0.037 7.012 
S10 Astragalus craibianus G. Simpson Legume P 0.025 6.854 0.012 —4.787 0.025 4.189 
S11 Tibetia himalaica (Baker) H. P. Tsui Legume P 0.015 —2.189 0.013 0.836 0.045 9.453 
812 Oxytropis ochrocephala Bunge Legume P - - - - 0.033 9.053 
813 Artemisia annua L. Forb A 0.104 7.983 0.008 —4.603 - - 
S14 Anaphalis lactea Maxim. Forb P 0.047 5.788 0.004 —8.326 0.045 8.129 
S15 Saussurea nigrescens Maxim. Forb P 0.026 6.589 0.014 1.703 0.062 10.708 
S16 Taraxacum lugubre Dahlst. Forb P 0.003 —10.802 - - 0.009 —5.799 
817 Cirsium periacanthaceum C. Shih Forb P 0.011  —10.812 - - - - 
Leontopodium nanum (Hook.f. & 
S18 Thomson ex Hook.f. & Thomson) Forb P - - 0.004 —8.314 - - 
Hand.-Mazz. 
gjo id à Pa da M Forb P P 2 : = 0.012 -5.540 
S20 Ajania tenuifolia Tzvelev Forb P - - - - 0.014 —0.940 
S21 Saussurea graminea Dunn Forb P - - - - 0.005 —9.958 
822 Potentilla anserina L. Forb P 0.033 7.290 0.011 —8.363 0.029 9.214 
S23 Potentilla nivea L. Forb P 0.004  —10.275 0.012 1.471 0.023 4.662 
s24. i tla cca Forb P — 0011  -8966 0.008  -8318 0.018 -5.468 
Cambess. 
S25 Trollius farreri Stapf Forb P 0.015 2.063 0.010 1.871 0.006 —9.927 
oe es me ia Forb P — 0.37 8548 0012 1.821 0.03 10286 
s27 demone — ad Forb P — 0019 1.721 0.005  -8314 0035 11.582 


To be continued 
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Continued 
No. Species Functional Life la 34a No mound 
group cycle — jv AOx — IV AO — IV AO; 
S28 . Gentianopsis paludosa (Hook. f.) Ma Forb A 0.007  —10.278 0.008 —4.527 0.007 —9.923 
S29 Halenia elliptica D. Don Forb A 0.029 6.296 0.005 -8.319 0.017 —1.326 
Sr Forb P — 0.008 -9.063 0007 4538 0.036 6.973 
Griseb. 
831 Gentiana leucomelaena Maxim. Forb A - - - - 0.014 2.500 
S32 Gentiana abaensis T. N. Ho Forb A - - - - 0.006 —10.050 
$33 Lancea tibetica Hook. f. & Thomson Forb P 0.011 —4.665 0.025 11.953 0.012 —0.878 
S34 Veronica eriogyne H. Winkl. Forb P 0.031 5.965 0.012 1.923 0.017 —5.860 
835 Pedicularis kansuensis Maxim. Forb A - - 0.004 —8.327 0.005 —9.936 
836 Epilobium royleanum Hausskn. Forb P 0.010 —9.068 - - - - 
837 Plantago depressa Willd. Forb A 0.020 2.053 0.024 5.604 0.010 5,725 
S38 Stellaria infracta Maxim. Forb P 0.027 8.462 0.037 14.024 - - 
839 Equisetum arvense L. Forb P 0.006 —7.658 0.006 4.123 - - 
S40 Geranium pylzowianum Maxim. Forb P 0.022 2.887 - - 0.009 —5.624 
S41 Allium sikkimense Baker Forb P - - 0.030 5.650 0.026 4.678 


Note: S1-S41 indicate plant species. - indicates no value. P, perennial; A, annual. 
mm | amound =33-4amound æ No mound 


(a) (b) (c) 


T-2JEOPOPOPJPJEEOE2E-30630 30303030 30:30 3030 3 ES 
NL ENO Io — UJ HMNDIBW ODN VAD INOS =VNWHAD ISOS 


(o 0/0 0/0 0/0 0,0 pd pepe pepe pepe pepe pe pepe pe pepe pepe pe pepe pe pepe pepe pe 008] 


Unununununununnun. 


0 1 2 3 4 5 6 700 04 08 12 16 20 24 2.80.0 0.4 0.8 12 16 2.0 24 28 
Shannon-Wiener niche breadth index Levins niche breadth index Hurlbert niche breadth index 


Fig. 2 Shannon-Wiener niche breadth index (a), Levins niche breadth index (b) and Hurlbert niche breadth 
index (c) of plant species under different types of zokor-made mound. S1-S41 indicate plant species. 


3.3 Niche overlap of plant species 


Most species pairs had niche overlap in 1 a (91.18%; Table S1), and 3-4 a mounds (87.19%; 
Table S2), and 100.0096 in no mound (Table S3). The proportions of species pairs with niche 
overlap of 70.6, 70.7, 70.8, 20.9 and 71.0 were 51.18%, 40.86%, 30.11%, 20.22% and 8.17% in 
1 a mound, respecitively; 47.41%, 40.52%, 36.21%, 32.76% and 19.33% in 3-4 a mound, 
respectively; and 55.35%, 46.97%, 38.59%, 27.45% and 9.18% in no-mound, respectively. The 
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total niche overlap in 1 a mound, 3-4 a mound and no mound was 58.32%, 59.31% and 65.91%, 
respectively. Niche overlap of plant species was initially small in new zokor-made mound and 
increased gradually as the succession progressed. The proportions of specie pairs in 3-4 a mound 
with niche overlap of >0.9 and >1.0 were 12.54% and 11.16%, which were 12.54% and 10.62%, 
and 5.31% and 10.15% higher than those of 1 a mound and no mound, respectively. These results 
indicate that competition among plant species increases over time, thus resulting in the increase of 
the interspecies niche overlap. 


3.4 Ecological response rate (R) of plant species 


According to the AO; we found 18 developmental and 13 declining populations in 1 a mound 
(Table 2). The main developmental populations were K. macrantha, E. nutans, R. tanguticus, S. 
infracta and A. annua and the main declining populations were C. periacanthaceum, A. fatua, T. 
lugubre, G. paludosa and P. nivea. In 3-4 a mound, we found 13 developmental and 16 declining 
populations. The main developmental populations were P. pratensis, S. infracta, E. nutans and L. 
tibetica and the main declining populations were P. anserina, P. kansuensis, A. lactea, H. 
elliptica, D. caeruleum, F. rubra, A. rivularis and L. nanum. In no mound, we found 17 
developmental populations and 17 declining populations. A. rivularis, S. nigrescens, R. tanguticus, 
T. himalaica, P. anserina, O. ochrocephala, F. rubra and A. lactea were the main developmental 
populations and F. ovina, G. abaensis, P. pratensis, E. nutans, S. graminea, P. kansuensis, T. 
farreri and G. paludosa were the main declining populations. 

The ecological response rates of different populations to the habitat are represented by the 
values of R. Among the developmental populations, T. distigmaticum (R-NBsw=2.41, 
R-NB,-0.86 and R-NBy=0.80), C. kansuensis (R-NBsw-1.13, R-NB,=0.40 and R-NBy=0.36) and 
A. rivularis (R-NBsw=1.04, R-NB,=0.38 and R-NBy=0.35) showed the higher R values, whereas, 
F. rubra (R-NBsw= —2.35, R-NB,- —0.81 and R-NBy= —0.73) was the most declining population 
with a larger R in 1 a mound (Fig. 3). In 3-4 a mound, T. himalaica (R-NBsw 71.60, R-NB; -0.49 
and R-NBy=0.40) was the main developmental population and C. kansuensis (R-NBsw= —1.97, 
R-NBi;- -0.63 and R-NBy= —0.56) and F. ovina (R-NBsw= —1.17, R-NB;- -0.38 and R-NBy= 
—0.33) (Fig. 3) were the declining populations. In no mound, T. distigmaticum was the main 
developmental population (R-NBsw=2.54, R-NB,-0.88 and R-NBy-0.81), and L. tibetica 
(R-NBsw- -1.82, R-NB;- -0.62 and R-NBy= -0.56), A. tenuifolia (R-NBsw = —1.71, R-NBi- 
—0.59 and R-NBy= -0.53) and H. elliptica (R-NBsw= —1.21, R-NBi7 —0.41 and R-NBy= -0.37) 
were the main declining populations (Fig. 3). 


3.5 Responses of niche breadth and R to zokor-made mound 


PCA ordination of niche breadth and R of plant species showed that the cumulative contribution 
rates of variance of the two principal components were 86.67% and 89.02%, respectively (Figs. 4 
and 5). Three types of mound were clearly separated by niche characteristics. Elymus nutans, P. 
pratensis, A. annua and S. infracta occupied larger areas in new mound (Fig. 4). E. nutans, P. 
pratensis, L. nanum, S. infracta and E. arvense appeared the most sensitive species in 3-4 a 
mound, three of which (E. nutans, P. pratensis and S. infracta) had a larger niche breadth 
compared with other species. In no mound, F rubra, K. capillifolia, O. ochrocephala, A. 
tenuifolia, S. graminea, L. stracheyi, T. lugubre, A. rivularis, G. leucomelaena, G. abaensis and G. 
ornata had a higher relevant with habitat compared with other species. Additionally, different 
plant functional groups such as FK rubra (grass), K. capillifolia (sedge), O. ochrocephala 
(herbaceous legume), A. rivularis (forb) and G. ornata (forb) all occupied a larger niche breadth, 
which indicated that the composition of plant functional groups was stable and the plant diversity 
was high in no mound treatment (Fig. 4; Table 2). 

PCA ordination of R in 1 a mound showed that E. nutans, P. pratensis, A. annua, A. rivularis, G. 
pylzowianum, P. depressa, T. farreri, H. elliptica and S. infracta were the most sensitive species 
after disturbance (Fig. 5). R values of A. rivularis, G pylzowianum, P. depressa and T. farreri 
were higher. In 3-4 a mound, S. nigrescens, P. nivea, R. tanguticus and V. eriogyne had higher R 
values, therefore they were the most sensitive species (Fig. 5). However, PCA ordination of R in 
no mound showed that K. macrantha, F. rubra, K. capillifolia, O. ochrocephala, A. sikkimense, G. 
leucomelaena and G. ornate had a higher correlation with habitat compared with other species 
(Fig. 5). Sedge K. capillifolia and forb G leucomelaena showed higher R values. 
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Fig. 3 Ecological response rates (R) of plant species under different types of zokor-made mound. R-NBsw, 
ecological response rate of Shannon-Wiener index (a); R-NB;, ecological response rate of Levins index (b); 
R-NBy, ecological response rate of Hurlbert index (c). S1-S41 indicate plant species. 
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Fig. 4 Principal components analysis (PCA) of niche breadth under different types of zokor-made mound. 
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Fig. 5 Principal components analysis (PCA) of ecological response rate under different types of zokor-made 
mound. S1—S41 indicate plant species. 


4 Discussion 


Disturbance by burrowing rodents is widespread across ecosystems. In tundra grasslands, rodents 
had positive effects on nutrient cycling by modifying the nutrient dynamics of plant communities. 
For example, disturbance increased nutrient content of forbs and grasses (Bon et al., 2020). In 
arid grasslands, plant community structure differed between mounds and landscape patches 
occupied by two rodents (Cynomys gunnisoni and Dipodomys spectabilis), and the rodents 
increased landscape heterogeneity, plant species richness and forb coverage (Davidson and 
Lightfoot, 2008). Similarly, disturbance by rodents (Rhombomys opimus, Cricetulus migratorius, 
Meriones meridianus and Dipus sagitta) altered the morphologies and nutrients of Haloxylon 
ammodendron in deserts, as well as increased soil organic matter, total nitrogen, available 
nitrogen and available potassium contents (Xiang et al., 2020). However, disturbance by rodents 
(Ctenomys mendocinus) reduced plant coverage, diversity and density (Lara et al., 2007). In 
tallgrass prairie, as selectively foraging on two legumes and one coneflower, voles altered plant 
community structure and reduced plant diversity (Howe and Brown, 1999). In grasslands, 
disturbance by rodents reduced productivity (Maron et al., 2014). Based on this study, we 
believed that disturbance by rodents also had a positive or negative effect on the plant community 
in alpine meadow. Plateau zokor, a typical alpine meadow rodent, plays an important role in 
ecosystem function, such as soil nutrient cycling, soil structure, vegetation composition and so on, 
creating heterogeneous habitat patches and being regarded as ecosystem engineers (Zhang and 
Liu, 2003). The mound-making by plateau zokor is an important determinant of plant community 
composition in the alpine meadow (Hu et al., 2017), however, how the zokor-made mound 
disturbance affected niche characteristics of plant community in alpine meadow had been often 
ignored. 

Niche breadth can be used as an indicator to measure the scale of species' utilization of 
environments or resources. Usually, the species with a higher importance value can use more 
resources or tolerate variable environmental conditions and distribute more widely, resulting in 
the larger niches (Dostál et al., 2016; Pannek et al., 2016). The larger the niche breadth of the 
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species, the more completely environments or resources are utilized, and the larger the species' 
habitat range (Sheth et al., 2020). Species differ in their competitive abilities to use habitat 
resources, affecting their niche breadths. In an ecosystem, the heterogeneous habitat created by 
the original community after disturbance causes the release of niche space. Then, the species 
affect community structure by competing or pre-empting niche space in heterogeneous habitat 
(Yang et al., 2018). Generally, because each species has an equal competitive opportunity to 
colonize bare zokor-made mound to occupy the niche during succession process, different types 
of forages could coexist, thus resulting in more abundant plant species in the climax community. 
Our study was similar to the results of Yang et al. (2018) and Niu et al. (2020), in that the 
composition characteristics of plant community were different at earlier stages of succession (1 a 
and 3-4 a) and inclined to become similar to surrounding undisturbed vegetation (no mound) as 
succession progressed. In no mound, the composition of functional groups of plant community 
was stable and plant diversity was high. Therefore, the later stage of succession in the alpine 
meadow after disturbance had higher species richness and plant diversity (Zhang et al., 2009; Niu 
et al., 2019). 

In this study, disturbance by rodents significantly reduced the species richness, coverage and 
Simpson diversity index (P«0.05; Table 1), but did not affect Shannon-Wiener diversity index and 
Pielou evenness index (P>0.05; Table 1). However, the height of plant community increased 
significantly (P«0.05; Table 1). The disturbance of zokor-made mound led to partial changes in 
the plant community structure, but there was no obvious damage to the overall structure of the 
plant community. The niche breadth of perennial grasses increased after disturbance (e.g., E. 
nutans and P. pratensis), and the response of plant community structure tended to increase the 
proportion of fine forages as succession progressed. The disturbance caused by zokor-made 
mound led to changes of soil nutrients and the competitive relationship between monocotyledons 
and dicotyledons plant communities, which promoted an increase of grasses (Zhang et al., 2004). 
The fibrous root structure of monocotyledons enables efficient nutrient absorption and utilization, 
thus making them more competitive on disturbed soil (Adler et al., 2007). Species differ in their 
ecological response to disturbance (Zhang and Liu, 2003; Wang et al., 2008; Niu et al., 2020), 
with perennial forbs usually responding aggressively (Zhang et al., 2003). In this study, K. 
capillifolia (sedge) and G. leucomelaena (forb) showed a stronger ecological response in 
undisturbed sites, which is consistent with the dominance of Kobresia meadow in the region 
(Wang et al., 2006; Li et al., 2011; Hu et al., 2017). It was similar to the result of Niu et al. (2020), 
in that the dominant species in the community were fine forages, and the proportion of sedges 
significantly increased at the end of succession. 

Analysis of competitive coexistence among species indicated that as succession progresses, the 
species that remained developmental included K. macrantha (perennial grass), S. nigrescens and 
R. tanguticus (perennial forbs; Fig. 6). These three species are the most common dominant 
species in alpine meadow communities of QTP (Li et al., 2011; Mu et al., 2015; Shi et al., 2015; 
Hu et al., 2017). In the present study, the niche breadth of these species was also higher than those 
of other species (Fig. 2). With the progress of succession, the species that transformed from 
declining to developmental were perennial grass F rubra, perennial sedge K. capillifolia, 
perennial legume T. himalaica and perennial forb P nivea and G. ornata (Fig. 6). The niche 
breadths of these species increased gradually with the succession (Fig. 2). Similar to our results, K. 
capillifolia and P. nivea were also reported to be the common dominant plant species with fine 
grazing capacity in the alpine meadow of QTP (Zhong et al., 2014; Guo et al., 2015; Xu et al., 
2018; Wang et al., 2019). The forb D. caeruleum and G. paludosa always remained declining with 
the progress of succession (Fig. 6), and their niche breadths were also smaller (Fig. 2). 
Additionally, perennial grasses E. nutans and P. pratensis, perennial forbs 7. farreri, V. eriogyne 
and G. pylzowianum and annual forbs H. elliptica and P depressa transformed from 
developmental to declining as succession progressed (Fig. 6), and niche breadths of these species 
decreased gradually over time (Fig. 2). This could be due to the interspecific competition 
interaction, which plays a vital role in determining community structure of either developmental 
or declining species, the interspecific competition became more obvious and the niche overlap 
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increased as succession progressed (Cutler, 2010; Wang et al, 2012; Yang et al., 2018). 
Interspecific competition became more intense, their niches tended to be more specialized and 
their niche breadths decreased. In the whole process of succession, most of annual plant species 
were in declining phase (e.g., A. fatua, A. annua, G. paludosa, H. elliptica, G leucomelaena, G. 
abaensis, P. kansuensis and P. depressa; Fig. 6). Annual species have poor competitive ability, 
and generally they are unable to compete against the more vigorous perennial species that have a 
stronger advantage in competition (Yang et al., 2018). Perennial monocotyledons, which have 
fibrous root system and strong ability to utilize soil nutrients, in particular have advantages in the 
heterogeneous environment following disturbance (Zhang et al., 2004; Adler et al., 2007). 
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Fig. 6 Coefficients of niche overlap of plant species under different types of zokor-made mound. (a), 1 a mound; 
(b), 3-4 a mound; (c), no mound. Bars are standard errors. * indicates significant difference among plant species 
at P<0.05 level. S1-S41 indicate plant species. 


5 Conclusions 


This study showed that niche breadth of different plant species diverged after disturbance by 
zokor-made mound. Mound construction didn't cause obvious damage to the overall structure and 
function of plant community in alpine meadow, but led to partial changes in plant community 
structure, which was conducive to stimulate the regeneration and vitality of plant community in 
alpine meadow. Mound increased the importance value and niche breadth of some perennial 
monocotyledonous forage species, which have more competitive advantages in the heterogeneous 
environment after disturbance. These perennial monocotyledonous species had obvious 
advantages in the utilization of heterogeneous habitat resources, and these perennial 
monocotyledons were often the main dominant species of alpine meadow. These dominant 
species with more extensive niche breadth had appropriate niche overlap with other species, 
ensuring their long-term ecological coexistence and adaptation. The composition of alpine 
meadow plant community would succeed in the direction of increasing proportion of fine forages, 
with the result that disturbance may have improved the grazing capacity of alpine meadow. 
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